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One of the outstanding problems in the numerical discretization of the Feynman- 
Kac formula calls for the design of arbitrary-order short-time approximations that are 
constructed in a stable way, yet only require knowledge of the potential function. In 
essence, the problem asks for the development of a functional analogue to the Gauss 
quadrature technique for one-dimensional functions. In PRE 69 (2004) 056701, it has 
been argued that the problem of designing an approximation of order v is equivalent to 
the problem of constructing discrete-time Gaussian processes that are supported on finite- 
dimensional probability spaces and match certain generalized moments of the Brownian 
motion. Since Gaussian processes are uniquely determined by their covariance matrix, 
it is tempting to reformulate the moment-matching problem in terms of the covariance 
matrix alone. Here, we show how this can be accomplished. 
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1 Introduction 

Since its introduction, the Feynman-Kac formula [1-3] has played a funda- 
mental role in the development of numerical algorithms capable of accounting for 
the physical properties of quantum systems made up of distinguishable particles. 
When utilized in tandem with the Monte Carlo integration technique, a powerful 
method is obtained: the path-integral Monte Carlo [4], which is capable of probing 
the quantum effects without any untestable approximations. The Feynman-Kac 
formula expresses the density matrix of a thermodynamic system as the expected 
value of a functional of the Brownian bridge 

p(a;,a;';/3) = ^^e-("'-"''/2"'Ecxp(-/3 / vlxriu) + aB°]du\ . (1) 

A second formulation of the Feynman-Kac formula is in terms of the full Brownian 
motion and reads 



{xle-f^^lip) = p{x,x';P)^{x')dx' = 

1 (2) 

= E exp(^-P J V{x + aBu)dujij{x + aBi) , 

where '0(x) is any square integrable function. In the above, p{x, x'; /3) is the density 
matrix for a one-dimensional canonical system characterized by the inverse tem- 
perature l3 — 1/(A:bT) and made up of identical particles of mass mo moving in the 
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potential V{x). The stochastic element that appears in Eq. {B^, < u < 1}, is 
a so-called standard Brownian bridge, defined as follows: if {Bu, tt > 0} is a stan- 
dard Brownian motion starting at zero, then the Brownian bridge is the stochastic 
process {B^, < u < = 0}, i.e., a Brownian motion conditioned on the 

event Bi = 0. As is well known, a Brownian bridge can be realized as the process 
— uBi, < M < 1} [5]. To complete the description of Eqs. (Q) and Q, we 
have Xr{u) — x + {x' — x)u and a = (^i^/J/toq)"'^^^- 

A problem of direct interest to the chemical physicist is the development of 
approximations supported on finite-dimensional probability spaces that have fast 
convergence for smooth enough potentials, as measured against the number of eval- 
uations of the potential function. Most desirably, such approximations should uti- 
lize only the potential function in their construction. For reasons of stability, they 
should converge, perhaps at a slower rate, for all continuous potentials that are 
bounded from below. Until recently, the fastest method available (as order of con- 
vergence) has been the trapezoidal Trotter discrete path integral method. The 
technique is usually derived by means of the Lie- Trotter product formula and 
an appropriate short-time high-temperature approximation. The formal asymp- 
totic convergence of the trapezoidal Trotter method and of related techniques was 
extensively studied by Suzuki [6, 7] for bounded operators and by Ichinose and 
Tamura [8, 9], among others. In particular, results of the last two authors [9] 
show that the symmetric Trotter-Suzuki approximation has optimal convergence 
0(l/n^) for sufficiently smooth potentials as far as pointwise convergence of their 
integral kernels is concerned. This type of convergence is also implied in the present 
paper. On the other hand, the non-existence theorem of Suzuki [6] makes it implau- 
sible that faster convergence can be achieved by utilizing short-time approximations 
constructed as functions of the kinetic and potential operators. 

Recently, a more general approach has been put forward by the present au- 
thor [10], who argued that, for sufficiently smooth potentials, there might exist 
direct short-time high-temperature approximations of arbitrary polynomial conver- 
gence order. Of course, these short-time approximations are generally not functions 
of the kinetic and potential operators. The construction of such approximations is 
based upon an "experimental" theorem on the pointwise convergence of the inte- 
gral kernels of the Lie-Trotter product formulas. Although not rigorously proved, 
this theorem seems quite plausible. The short-time approximations considered are 
based on carefully designed finite-dimensional approximations to the Brownian mo- 
tion entering the Feynman-Kac formula. Basically, the Brownian motion is replaced 
by some discrete-time Gaussian process that is supported on a finite-dimensional 
probability space. A set of functional equations involving some generalized mo- 
ments of the Gaussian process have been shown to control the order of convergence 
u. Because the number of equations increases in an exponential fashion with ly, 
explicit solutions have been obtained only for v = 4. 

In the present work, we exploit the fact that both the Brownian motion and 
its replacement are Gaussian processes and are, therefore, uniquely determined by 
their covariance matrices. We thus show how to express the functional equations 
in terms of the covariance matrices alone. Hopefully, the new equations will prove 
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easier to utilize in a complete mathematical proof of the existence of short-time 
approximations of arbitrary order. In the appendix, we give a general convergence 
theorem regarding the construction of finite-dimensional discrete approximations 
to the Feynman-Kac formula. 



2 Statement of the moment-matching problem 

Perhaps one of the oldest techniques for simulating a Brownian bridge (or Brow- 
nian motion) is via random series. As such, let {Xk {T)}k>o be any orthonormal basis 
in L2[0, 1] such that Ao(t) ee 1, let 



Afe(u) 



Afe(T)dr for A: > 



and let a :— {oq, ai, . . .} be a sequence of independent identically distributed stan- 
dard normal random variables. By the Ito-Nisio theorem [11], the random series 
Sfe>o ^*:Afe(u) is uniformly convergent almost surely and equal in distribution with 
a standard Brownian motion starting at zero. By the construction of a Brow- 
nian bridge as the process — uBi, < it < 1} and the fact that Ao(u) = u, it 
follows that X]fc>i O'feAfe(u) = X)/c>o ~ '^o^ equal in distribution with a 
standard Brownian bridge. If fl is the set of all sequences a := {qq, oi, . . .} and if 

oo 

dP[a] = Y[ dfi{ak) with dfi{z) = {2n)-'^/^ exp(-2V2)dz 

k=0 

is the probability measure on Q associated with the sequence of independent random 
variables a := {oq, ai, . . .}, then the Feynman-Kac formula given by Eq. reads 



V 



k=l 



du 



(3) 

Here, pfp{x,x'; P) = exp[— (x' — x)^ /2a^]/ {2na^y^^ is recognized as the density 
matrix of a free particle. The alternative formulation given by Eq. ^ becomes 



p{x,x';P)= / dP[a]cxp\-f3 V x + cr^^ akAkiu) dul 'ip{x + aao) ■ (4) 



k=0 



Eqs. (O and Q are appropriately called random series representations of the 
Feynman-Kac formula in the chemical-physics literature [12, 13]. 

Eq. ||2J) is suggestive of some sort of numerical approximation to the Feynman- 
Kac formula, namely 



Pn{x,x';P)=pfp{x,x';f3) j dP[a]exp <^ -(3J2wiV Xr{9i) + (jJ2<^kJ^k{di) 

i=l k=l 



(5) 
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Here, the non-negative weights Wi (we assume Wi = 1) and the knots 6i define 
a quadrature technique on the interval [0,1]. It goes without saying that such 
approximations are convergent as — > c» and —^ oo under mild assumptions 
for the potential V{x): boundness from below and continuity. In order not to 
disrupt from the flow of the presentation, we give the simple proof in the appendix 
(see Corollary Thus, we are in no shortage of quadrature formulas. We just 
want the faster ones. 

Another utilization of Eq. Q is as a short-time approximation in a Lie-Trotter 
product of the form 

Pn{x,x']f3)= I dxi... [ dxnpi''Ux,xi;^^]...plj''Uxn,x';—^]. (6) 

The gain in interpreting Eq. Q as a short-time approximation is that the re- 
quirement that the functions {Afc(u); 1 < fc < n^} are constructed according to the 
Ito-Nisio prescription can be relaxed. We denote this by utilizing a tilde overscript, 
so that the short-time approximation reads 



p'l)''\x,x';(3) = pfp{x,x';P) I dp{ai) . . . j d^(a„Jx 



X exp < —(3 / V 







Xr{u) + afeAfc(w) 



k=l 



dLu{u] 



(7) 



The functions {Afe(u); 1 < k < n^} are required to be continuous (that is, bounded 
on the quadrature knots) and vanish at both end points. At this moment, we 
should emphasize that this special choice of short-time approximation is made 
while bearing in mind its usefulness in Monte Carlo simulations. Due to some 
special properties, especially the availability of the fast sampling algorithm [14], 
the subsequence of Lie-Trotter products with n = 2*"' — 1 is of utmost practical 
interest. A simple proof of the convergence of this Lie-Trotter subsequence is given 
in the appendix (see Corollary , again for continuous and bounded from below 
potentials. 

A second requirement that we ask of the construction given by Eq. (Q has 
to do with the symmetry of the density matrix, which should reflect itself in the 
symmetry of the short-time approximation. As such, we require that the discrete 
probability measure 

duj{u) — Wid{u — 9i)du (8) 

i=l 

defining the quadrature technique on [0, 1] must be symmetric about i. Also, the fi- 
nite dimensional process X]fe=i Ofe-^fe(u) must be invariant under the transformation 
u' — 1 — u. That is, we require the equality in distribution 

^afeAfe(u) = ^afeAfe(l - u) . (9) 

k=l k=l 
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The time symmetry of the process be enforced, for example, 

by restricting the functions {Afc(u); 1 < A: < n^} to the class of symmetric and 
antisymmetric functions. To understand these requirements, notice that the Her- 
miticity of the density matrix stems from the symmetry of the Lebesgue measure 
on [0, 1] as well as from the time symmetry of the standard Brownian bridge S", 
i.e., the fact that < u < 1} is a Brownian bridge equal in distribution 

to {5°, < u < 1}. Because the random sum J2t=i'^k-^k{u) is intended as a 
replacement for the Brownian bridge, it is convenient to introduce the notations 



n 

Bl = ^akAk{u) and 13^ = Qqu + for < m < 1 . (10) 

fc=i 



Thus, the Brownian bridge (motion) is approximated by a simple Gaussian process 
that is supported on a finite-dimensional probability space. Moreover, only the 
values of the process for the discrete times represented by the quadrature knots 9i , 
02, ■ • • , drig are relevant for the construction of the short-time approximation. 

The symbol (v) appearing as a superscript in Eq. 10 denotes the order of 
convergence of the short-time approximation. This is defined as the largest v for 
which the so-called convergence operator T^^ expressed by 

/m Po'^'la^-a:';/?) - p(x,x';/3) ip{x')dx' 



(T.^)(.) ^hm ^ . (11) 

is well defined at least on the class of infinitely differentiable and compactly sup- 
ported functions ip{x). In Ref. [10], it is claimed but not rigorously demonstrated 
that 

lim {n + iy[p,,{x,x';P)-p{x,x';/3)]^P''+^ C (x c~spHj.^^-{i-e)pH ^ 

n^oc Jq \ I 

(12) 

where Pn{x,x'; /3) is defined by Eq. A rigorous proof of Eq. (|12() is beyond 
the mathematical abilities of the present author. Based on the more or less formal 
arguments presented in the aforementioned reference, it is very likely that the 
statement is true. The author would be very grateful to the mathematically more 
inclined reader who may want to investigate the problem and prove or disprove the 
assertion (again, the case rt = 2*^ — 1 suffices for all practical purposes). 

Nonetheless, Eq. H12(l states that the convergence of the Lie-Trotter product 
is as fast as l/n", which explains the nomenclature regarding the order of conver- 
gence. We should emphasize that merely the existence of the convergence operator 
expressed by Eq. Hll() sets some constraints on the smoothness of the potential 
function. The natural class of potentials to study the problem of constructing 
short-time approximations of arbitrary order is the class of continuously and in- 
finitely differentiable functions for which 



1 



e-^V(2a) 



dz < oo , (13) 
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for all X e M and a > and for all integers fc > and j > 1. This condition is 
necessary in order to ensure that we recover the original potential V{x), derivatives 
V^'^^x), or products of such functions from their Gaussian transforms, in the limit 
that a — > (see Theorem 3 of Ref. [15]). We mention that this class of potentials 
does not include some pathological infinitely differentiable and bounded from below 
potentials such as exp(x'*) or even cos[exp(a;'')]. 

In these conditions, according to Theorem 4 of Ref. [10], the convergence oper- 
ator exists if and only if 

E [{B.y'iMoy^M.y' . . . (M2^_2)-''^-] = e [{Biy^iMoYHMiy' . . . (a12^_2)^'^"" 

(14) 

for all 2/i-tuples of non-negative integers (ji, j2, . • ■ , j2/i) such that X]fe=i ^Jfe — 2/i 
and 1 < n < V- The random variables Mk and Mk are defined by 



M, 



k 



/ (B„)'=du and Mk= (B„)'=dc^(u) = V u;,(BeJ^ (15) 
Jo Jo 



respectively. 



3 Reformulation of the generalized moment conditions in terms of 

covariance matrices 

Since they are Gaussian processes with continuous paths, the Brownian motion 
and its approximation are uniquely determined by the covariance matrices 

oo 

7(u,t) =E(B„S,) =ur + ^Afe(w)Afe(r) =min{u,T} (16) 

and 

7(u, r) = E{BM ^UT + J2 h{u)Ak{T) , (17) 

fc=i 

respectively. As such, at least in principle, the relations given by Eq. (|14|l can be 
formulated in terms of these covariance matrices alone. In this section, we show 
how this can be done. 

Let J2f_i, denote the set of solutions of the Diophantine equation X]fe=i ^Jk = 2/i. 
For each ( G J2m, we define the integer 

d(C) =j3+j4 + ---+j2M (18) 

and 

"(C) = [ji + J3 + 2j4 + • • • + (2Ai - 2)j2m] = - J2 - d{C) , (19) 

respectively. We also define the differentiation functional I?^ acting on the space of 
infinitely differentiable functions /(Aq, Ai, . . . , ^d{c)) that associates to each / the 
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following partial derivative evaluated at the origin 



— /(0,0,...,0). 



dX 



2/^-2 

rf(C)-i2„ + l'^^"'i(C)-j2^,+2 



dxl'f; 



■ dx 



2/^-2 



The differential order of is ji + + 2j4 + • • • + (2/i — 2)^2/^ = 2n(C). 

The reader will understand the need for this rather cumbersome notation short- 
ly. Going back to Eq. let us notice that the equality 



Idu = 1 



lAuj{u) ^ Mq 



means that the factors containing AIq and Mq cancel out. Next, we utilize the 
identities 



A=0 



and Mk — kl I du-—r e 
dA*^ 



A=0 



to express the first factor in Eq. (|14(l as 



ji!(l!p...((2/.-2)!) 



V dA-Ji 



A=0 



du 



dA2M- 



A=0 



A=0 



(21) 



Expanding the parenthesis and interchanging the order of integration and differen- 
tiation, we obtain the result 



where 



/c(^o, Ai, . . . , Xd{Q) 



/■I /■! /d(C) \ 

f dui--- dud(c)Eexp J2 ^i^u, ] 

Jo \ 1=0 / 



and Mo = 1. 

To evaluate the function f(^, one may utilize any random series and compute 



Eexp( J2^iBu. =Eexp J2 Y.<^k^k{ui) 
, i=o / V 1=0 fc=o 



fc=0 



XiKk{ui) 

1=0 
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The last term equals 



1^(0^(0 



2 i=o j=o 



exp <J E E ^^^j { E Afc(M,)Afc(uj) ) [> = exp <{ - E A,Aj7(u,, u^) 

, fc=0 



irf(C) <i(C) 



2 i=o j=o 



Therefore, the explicit expression of the function is 



Jo I ^ i=0 j=0 J 



/c(Ao, Ai, . . . , Ad(^)) 



We remind the reader that uq = 1. 

The function can be replaced by some (1 + rf(C))-dimensional polynomial of 
degree 2n{C,). Indeed, starting from the expansion 



<i(c)<i(C) 

E E AiAj7(ui,Uj) 

j=0 j=0 



one observes that, upon the action of D^, the differential order of which is 2n{C,), 
only the term with k = 't-(C) survives. The terms of lower degree are washed out 
by the process of differentiation, whereas the terms of higher degrees cancel when 
the respective derivatives are evaluated at the origin. Thus, the left-hand side of 
Eq. (|14|l is given by the expression 



ji!(l!r^---((2/i-2)!)-'"^- 
2«(C)n(C)! 



Dcfc, 



where 

/c(Ao, Ai, . . . , \d{o) 



dui ■ 



du 



d(C)d(C) 

E E AiAj7(w«,Wj) 

1=0 j=0 



AO 



(22) 



(23) 



In an analogue manner, one demonstrates that the right-hand side of Eq. (|14() 
can be written as 

ji!(l!)^3...((2^_2)!)-'--- 



2"(C)n(C)! 



Dcfc , 



(24) 



where 



/c(Ao,Ai,...,A<j(c)) = / duj{ui)--- dLu{ud(Q) 



d(Od{0 

E E A»Aj7(ui,Uj; 

i=0 j=0 



<0 



(25) 

Comparing Eqs. (|22|l and 124|l . we see that Eq. (|14|) is equivalent to the equality 

DJc = Dcfc (26) 
with and defined by Eqs. (^5)1 and respectively. 
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As Eqs. (I22|l and (|24|l show, the polynomials ft; and /(; depend only on some 
integers d{() and n{C) that have the property 

d{0+n{0 <ti<iy. 

Thus, a sufficient condition that the equality expressed by Eq. (|2()|l holds for all 
C G J^fj. and 1 < fi < v is that the equality of polynomials 

fn,d{\), Ai, . . . , \d) 
= / dui 

"'0 

holds for all integer d and n such that d + n < v. If c?i < o?2 then the polynomial 
of dimension 1 + di is a particular case of the polynomial of dimension 1 + d2, with 
the last ^2 — di lambda's set to zero. Thus, it is enough to check the above equality 
for the cases with d + n = v. 

Let us show that checking the equality expressed by Eq. H27|l for all integers n 
and d such that n + d — v is, also a necessary condition. Notice that the polynomials 
appearing in Eq. H27|l are symmetric under the permutation of the variables Ai, 
A2, . . . , \d- This is due to the symmetry of the covariancc matrices 7(1*, t) and 
7(u,r). In these conditions, the differentiation functional given by Eq. (|20(l is 
the most general differential expression of order 2n, provided that we let the index 
C = {jii j2, ja, ^4, • • ■ ,j2v\ lie in the set S that contains all possible indexes for which 

j2 = 0, j3 + j4 H h j2v = d, and ji + js + 2j4 H h (2i/ - 2)^2,. = 2n. Following 

some permutation of the variables Ai, A2, . . . , A^, any other differential expression 
can be obtain from and, by the aforementioned symmetry of polynomials, is equal 
to a differential expression of the type given by Eq. H2U|) for some j E S. Therefore, 
the equality V^fn d — T^cfn.d for all C G 5* implies fn,d — fn.d- It remains to prove 
that, as C spans J2,y, it also spans S. Clearly, any index C & S has the property 
ji + 2j2 + 3^3 + • • • + 2iyj2i, — 2(n + d) — 2z/. Therefore, S C ^21/ and the claim of 
necessity is proved. 

We summarize the results of the paper in the following proposition. 

Proposition 1 A short-time approximation of the type given by Eq. ^ has con- 
vergence order v if and only if the equality of polynomials expressed by Eq. \21^ 
holds for all integers d and n such that d + n = v. 

The above proposition can be utilized together with the multinomial formula 
to generate the necessary conditions that the trial covariancc matrix must satisfy 
in order for a short-time approximation to have convergence order v. Most likely, 
more useful statements can be obtained by the mathematically more inclined reader. 
Nevertheless, we have achieved our goal of formulating the conditions in terms of 
the covariancc matrices alone. 

This work was supported in part by the National Science Foundation Grant No. CHE- 
0345280 and by the Director, Office of Science, Office of Basic Energy Sciences, Chemical 
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1 = 7 = 



d d 

E EAiAj7(M^,UjO 
i=0 j=0 
-1 n 



(27) 
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Appendix 

In this appendix, we establish an easily provable theorem that is useful in many 
situations and provides the basic mathematical background for the discretization 
of the Feynman-Kac formula. We then utilize the theorem to demonstrate several 
statements made in the text. 

To begin with, let us formalize what we mean by quadrature rules. By a quadra- 
ture scheme on [0, 1], we understand a sequence of pairs of vectors of variable lengths 
g„, namely the weights {wn.ijWn.2i ■ ■ ■ ,Wn^q„} and the knots < u„_i < u„^2 < 
• ■ • < Un.q„ < 1, constructed such that the Wn/s are non-negative for all 1 < i < qn 
and n > 1 and such that 

lim Wn,ih{un,i) ^ / h{u)du (28) 

n — »oo ^ ^ I r, 

1=1 •^'J 

for all continuous functions h : [0, 1] —>■ M. The mid-point, trapezoidal, Simpson, 
and Gauss-Legendre rules are well-known examples of quadrature schemes. The 
non-negativity of the weights is a stability requirement that simply says that it is 
not all right to get a negative answer if a positive function is integrated. 

Let {5„(w; u); < u < l}„>i be a family of random processes with continuous 
paths supported on a probability space (f2, J-", P) and such that 5„(w; u) converges 
to some limit i?2(w) uniformly in u, P-almost surely. Assume that the limit B^(uj) 
is equal in distribution to a Brownian bridge. Given an arbitrary quadrature rule, 
we can define a sequence of approximations to the density matrix by the prescription 



Pn{x,x';P) ^ pfp{x,x'; f3)Ee^pi^-PJ2'^riAV [xr{un,i) + aSniuj;un,i)] |, (29) 

for n > 1. We shall call the above prescription the standard discretization of the 
Feynman-Kac formula. The nomenclature is motivated by the following funda- 
mental theorem. 

Theorem 1 // V{x) is continuous and bounded from below, then the sequence 
Pn{x,x'; P) is bounded by C{l3)pfp{x,x'; (3), for some positive C{j3) < oo, and 

Pn{x,x'; P) ^ p{x,x'; P) , as n ^ oo . (30) 

Proof. Let Vq be a non-positive lower bound for V{x). Eq. H28|l specialized for 
h{x) = 1 says that the sequence X]i=i '^^n,i is convergent, thus bounded, say by c > 
0. Then the integrand of Eq. (|29|) is bounded from above by C{/3) — exp(— c/?Vo) 
and we have 

Puix, x'; P) < C{P)pfp{x, x'- P) . (31) 
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The remainder of the theorem follows from the dominated convergence theorem 
and the Feynman-Kac formula, as soon as we prove that 



lim Wn,iV [Xr{Uni) + CrSniuj; Un,i)] = / V[xr{u) + aB° 



(w)]du, a.s. (32) 



As such, let w G be given. Since x.r(u) + oB^iio) is almost surely contin- 
uous in M — as uniform limit of continuous functions Xr{u) + aS„{Lo;u) — it is 
bounded by some constant M > 0. In fact, by uniform convergence, we have 
\xr{u) + aSnito;u)\ < 2M for large enough n. Since V{x) is continuous, it is uni- 
formly continuous on the compact set < 2M . Consequently, given e > 0, there 
is ry > such that |y(a;) — t^(y)| < e whenever — j/| < rj. Nevertheless, by uniform 
convergence, 

I [xr{u) + aSn{uj; u)] — [xr{u) + cri?"(cj)] | < ?7 
for large enough n and so, 

\V[xriu) + aSniio\u)] - V [xr{u) + aB^{Lj)] I < e. 

Since J2i=i '^n,i < c, it follows that 

?n In 

^Wn,iV [Xr{Un,i) + <TSn{uJ; Und)] ~~'^Wn,iV XriUnd) + Cri?.°,^^ (oj) 



< ce. 



(33) 

By Eq. H28|) and the continuity of V \xr[u) + (jB^[uj)\ as a function in u, we also 
have 



qn 



V[Xr{u) + aB^{uj)]du- '^Wn^iV Xr{Un,i) + fji?"^^ (w) 



for n sufhciently large. Combining Eqs. and we obtain 

'S^Wn,iV[Xr{Un,i) + (TSn{i^]Un,i)]- / V[Xr{u) + l7Bl{uj)\du 
.-1 "'O 



< e 



(34) 



< {l + c)e. 



Since e is arbitrary, the almost sure convergence appearing in Eq. (|32|) is demon- 
strated and the proof of the theorem is concluded. □ 
Observation. Because pn{x,x';(3) is bounded by C{(3)pfp{x,x'](3), the dominated 
convergence theorem and the above theorem also imply convergence in the strong 
topology, that is, 

Pnix^x' ; (})^j{x')dx' / p{x, x' ; (3)ip{x')dx' 



for all square integrable ^{x). In fact, by choosing the sequence of processes 
Sn{u:,n) to be constant and equal in distribution to S", Th.^ produces various 
versions of the Trotter convergence theorem for various quadrature rules. 

As everywhere else in this paper, in the following, it is understood that the 
potential V{x) is continuous and bounded from below. 
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Corollary 1 // pn{x, x' \ /3) is defined by Eq. 0), then pn{x, x'; 13) — > p{x, x'; 13), as 
n ^ oo. 

Proof. Follows from Th. ^ and the uniform convergence of the random series 
Sfc>i '^k^kiu) to a Brownian bridge, as guaranteed by the Ito-Nisio theorem. □ 

Corollary 2 // pn{x, x' \ (3) is defined by Eq. then p2k_i{x, x' ; j3) — > p{x, x' ; (3), 
as fc —> oo. 

Observation. The proof we construct for the subsequence n = 2'^ — 1 is based on 
a special form of the Lie-Trotter product that is called the Levy-Ciesielski form. 
The information from the paragraphs below is taken from Ref. [10]. In fact, for 
n = 2^^ — 1, any Lie-Trotter product can be put in a Levy-Ciesielski form, as shown 
in Ref. [14]. This form has important advantages when it comes to the Monte 
Carlo implementation, such as fast computation of paths and fast sampling. For 
this reason, we can restrict our attention to the subsequence n = 2*^ — 1, without 
any loss of generality for actual applications. 

Proof of the corollary. Let {a/ 1 < Z < fc, 1 < j < 2'~^} and {bij; 1 < I < 
g, 1 < j < 2*^} be two independent sets of i.i.d. standard normal variables. Let 
{Fi j{u)\ Z > 1, 1 < j < 2'~^} be the system of Schauder functions on the interval 
[0, 1]. The Schauder functions are the primitives of the Haar i^([0, 1]) wavelet bases 
and can be generated by contractions and translations as follows. Let Fi i{u) : M 
M be defined by 

r u, 7. e [0,1/2], 

Fi,i{u) =1 l-u, [1/2,1], (35) 

[ , elsewhere . 

Then, 

F;,,H = 2-('-i)/2^^i,i(2'-iu-j + l) (36) 

for alH > 1 and 1 < j < 2'~^. Extend the functions {A/(m); 1 < Z < n^} outside 
the interval [0, 1] by setting them to zero [the same way the first Schauder function 
Fi.i(m) was extended to the whole real axis in Eq. 1)35(1 ] and define 

Gzj(u) ==2-'^/2^,(2'=u-j + l) (37) 

for 1 < Z < rii, and 1 < j < 2^. 

Let du}k{u) denote the discrete measure associated with the quadrature scheme 
specified by the nq2^ (not necessarily different) quadrature knots 

u[^^^2-\e,+j~\), l<i<nq, 1<J<2'= (38) 

and the corresponding weights 

w[^j = 2-''w, . (39) 

The new quadrature knots u'^ j are obtained by contractions and translations of the 
original knots 6i . 
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With the convention that a; 2i-i+i — and 6;. 2^+1 — for aU I £ l,k, we have 



Pnix,x';P) 



fc 2' 



X y"d6i,i...y"d6„„,2'= (2^)"*"+'^"''/'exp(^ 



1 



exp -^E E^^L- 



X exp< — /? / V 



Xr[u) + crEa;,[2'-i«]+i -P/^p'-iuJ+iCw) 



1=1 



dwfc(u) 



(40) 



where [2'~^m] and p'^u] are the integer parts of 2'~^u and 2'^m, respectively. 

Let us verify that the sequence of discrete measures dujk (u) defines a quadrature 
scheme. If h : [0, 1] ^ M is a continuous functions then, by its uniform continuity, 
for any e > 0, there is > such that \h{u) — /i(t)| < e whenever |m — r| < 77. 
Now, pick k large enough that 1/2*'' < rj. The Lebesgue integral of h over [0, 1] can 
be broken in 2*"' smaller parts over the intervals [{j — 1)72*^, j/2'^], for 1 < j < 2*^. 
Since the length of these intervals is smaller than rj, we have both 



/(j-l)/2'= 



h{u)du-2-''h{j/2'') 



< 2-'=e 



and 



J^wl^^h{u'J-2-''hij/2') 



i=l 



< 2-^e. 



(41) 



(42) 



The latter inequality is true because the points u'^ j are in the interval [(j 
l)/2'', j72*''] and because 



i=l 



From Eqs. ||1T|| and we learn that 



/ /i(u)du - Vt«^ /i (u^ 

J(7-i)/2'= TTT 



and, by summing over all j, we obtain 

«/2' 



< 2 • 2-^t 



(43) 



2* 

E 



/ /i(u)dM-Vw^ 

iO-l)/2'= ^ 



< 2e. 



(44) 
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The left-hand side of Eq. I|44|l is clearly larger than 



E 



/ /i(u)d-u - Vw- « ) 



h{u)du — / /i(u)dci;/c(u) 



Since e is arbitrary, it follows that 



1 



lim / h{u)du!k{u) =^ / h{u)du 
k^^Jo Jo 

and since h is arbitrary, it follows that the sequence of discrete measures dujk(u) 
defines a quadrature scheme. 

The first part of the random series appearing in Eq. (|40() is the Levy-Ciesielski 
series, which converges uniformly to a Brownian bridge. The corollary readily 
follows from the result in the preceding paragraph and Th. ^ if we prove that the 
tail series 

Tk(b;u) = ^6;,[2'=„]+i G;j2'»n]+i(w) (45) 
1=1 

converges to zero uniformly almost surely, as fc — > oo (in other words, if our cor- 
rection term does not ruin the uniform convergence of the Levy-Ciesielski series). 
Let M > be a common bound for the functions {Afe(u); 1 < A; < rii,}. Then, with 
the help of Eq. H37|l , we compute 



Taking the expected value, we get 



3 

1=1 j=l 



E 



max|T,(M)| <^(3n.2^) = ^ 



Now, Chebyshev's inequality produces 

max |T,(M)| > e) < ( max |T.(M) |) ' < ^ 



and so. 



^ max \nib; -)| > . j < ^ E ^ = ^ < - ■ 

k=l ^ ^ ^ ' k=l 



The first Borel-Cantelli lemma implies that 



Pi max \Tk(b;u)\ > e i.o. =0, 

* 0<M<1 I ^ ' 'I - ' 
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which means that, with probabihty one, there is a rank K >1 such that 

max Tfc (6; m) < e , M k> K . 

0<-ii<l ' ' 

Letting e go to zero through the countable sequence €j = we obtain the almost 

sure uniform convergence to zero of the tail series T/j(6; u). The proof of the corollary 

is concluded. □ 
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